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ABSTRACT 

The nature and even the existence of a putative planet-mass companion (" Fo- 
malhaut b") to Fomalhaut has been debated since 2008. In the present paper 
we reanalyze the multi-epoch Hubble Space Telescope (HST) optical images on 
which the discovery claim was based. We confirm that the HST images do reveal 
an object in orbit around Fomalhaut but the detailed results from our analysis 
differ in some ways from previous discussions. In particular, we do not confirm 
flux variability over a two-year interval at 0.6 /im wavelength, we detect Fomal- 
haut b for the first time at the short wavelength of 0.43 /im, we find that the HST 
image of Fomalhaut b at 0.8 /im may be extended beyond the PSF, and we can- 
not determine from our astrometry if Fomalhaut b will cross or not the dust ring. 
The optical through mid-infrared spectral energy distribution (SED) of Fomal- 
haut b cannot be explained as due to direct or scattered radiation from a massive 
planet. We consider two models to explain the SED: (1) a large circumplanetary 
disk around a massive, but unseen, planet and (2) the aftermath of a collision 
during the past 100 years of two Kuiper Belt-like objects of radii ~ 50 km. 

Subject headings: Methods: data analysis, Methods: observational, Tech- 
niques: high angular resolution, Techniques: image processing, : planetary 
systems. 
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1. Introduction 

Direct imaging is the appropriate tech- 
nique for the study of exoplanets with 



nomical units ( 


Marois et al. 




200? 




2010 


Kalas et al. 




2008: 


Laeranee et al. 




20091) 



As the planetary atmospheric thermal 
emission or scattered light is detected, 
detailed multi-band photometry or spec- 
trometry can be used to probe the at- 
mospheric composition and physical prop- 
erties. However, these studies are chal- 
lenging given the very high contrast and 
small angular separation between a star 
and planet. In some systems the presence 
of a planet before it is detected directly 
can be suggested by the geometry of a 
circumstellar debris dis k. For e x ample , 



Wvatt et all (ll999l). iKalas et al.l (j2005h 



and lQuillen et al.l ( 120061 ) had predicted the 
likely exi stence of a planet aro und Fomal- 
haut and iMouillet et al.l ( 119971 ) of a planet 
around /3-Pictoris. 

In t he case of Fomalhaut (440±40Myr , 
7.7 pc. lMamajekll2012l : [Van Leeuwenll2007l ) . 
a candidate plan et was announced by 



Kalas et al.l (120081 hereafter K08). Sur- 
prisingly, the candidate was not detected 
in deep near infrared images in H and Lp 
bands, but rather in Hubble Space Tele- 
scope (HST) images in the visible where 
planets are not expected to emit much 
thermal light. The K08 planet model that 
best fit the 2008 photometry is a < 3Mj 
Jovian planet surrounded by a large cir- 
cumplanetary disk; the observed optical 
light is mostly scattered by the disk, not 
by the planet itself. In this model, it is as- 
sumed that Ha emission (dust accretion or 
hot planetary chromosphere) explains the 
unusual 50% variability of the reported 



flux at 0.6 /im over a two year time inter- 
val. B ased on the i r astr ometric mesure- 
ments, IKalas et al.l ( 120101 ) also announced 



that Fomalhaut b is likely to cross the dust 
ring. 

A few years later, as they could not 
detect the object a t 4.5 /im with Spitzer, 



Janson et al.l (120121 ) concluded that "there 
is almost certainly no direct flux from a 
planet contributing to the visible-light sig- 
nature" and they proposed an optically 
thin dust cloud with or without a central 
object in the super-Earth regime to explain 



the K 08 photometry. Kennedy &: Wvatt 
( 1201 ll ) also rejected direct detection of 
massive planets and explain the photome- 
try at 0.6-0.8 /im to be a consequence of a 
swarm of satellites around a 2-100 ME ar th 
planet. 

Motivated by the controversial status 
of Fomalhaut b within the community, 
including even doubts of its actual exis- 
tence, we decided to conduct an indepen- 
dent analysis of the HST public data that 
were recorded in 2004, 2006 and 2010. Af- 
ter describing the observational method 
in §[2] and our data reduction in §E1 we 
analyze the images to confirm that Foma- 
lhaut b is a convincingly real detection 
and that it is gravitationally bound to the 
star. In § 14.21 we study various possible or- 
bits to determine if the current astrometry 
can confirm or reject a dust belt cross- 
ing trajectory (su ch as one announced by 
Kalas et al.ll2010l ). In §H]we estimate the 



object's photometry and possible origin as 
a circumplanetary disk around a planet or 
the aftermath of a collision of two Kuiper 
Belt-like objects, while considering that 
Fomalhaut b is not ( § 14 . 3 j) or is (§ 14. 4p spa- 
tially resolved in the HST images. 
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ACS/HRC 


1.8 


F435W 


9 


6525 


3 


5.8 


2006 Jul. 15-16 


ACS/HRC 


3.0 


F435W 


9 


6435 


3 


5.8 


2006 Jul. 19-20 


ACS/HRC 


3.0 


F606W 


28 


7240 


4 


6.0 


2006 Jul. 18 


ACS/HRC 


1.8 


F814W 


20 


5280 


3 


6.0 


2006 Jul. 19 


ACS/HRC 


3.0 


F814W 


27 


4942 


3 


6.0 


2010 Jun. 14 
2010 Sep. 13 


STIS/50CORON 
STIS/50CORON 


2.5 
2.5 


CLEAR 
CLEAR 


3 
16 


630 
3000 


7 


157.0 



Table 1: Fomalhaut observing log. Column "Im." gives the number of useful images. Column 
"Exp." is the total integrated time of the useful images. "Roll" is the number of roll angles 
in the sequence and "FOV rot." gives the total FOV rotation during the sequence. 



2. Observations 

The data that we consider in this pa- 
per were obtained with HST with the Ad- 
vanced Camera for Surveys (ACS) in 2004 
and 2006 (programs 10390 and 10598), 
and the Space Telescope Imaging Spec- 
trograph instrument (STIS) in 2010 (pro- 
gram 11818). The ACS data were ac- 
quired with the High Resolution Chan- 
nel (HRC) in its coronagraphic mode 
with 1.8" and 3.0" focal plane occult- 
ing masks and the F435W (430 ± 50 nm), 
F606W (595±115nm), and F814W (825± 
115 nm) filters. For the STIS data, the 
50CORON configuration was used with its 
clear aperture (600 ± 220 nm). For all se- 
quences, images at several roll angles were 
taken so that the stellar diffraction pattern 
can be subtracted while keeping the flux of 
any point sources. Tab. [1] gives the dates 
of the observations, the instrument config- 
urations (filter, coronagraph) , the number 
of useful images with the corresponding 
integration time, and the number of roll 
angles as well as the total rotation of the 



field of view. 

3. Data reduction 

3.1. ACS 

We start from the drz drizzled images 
produced by the ACS pipeline (geometric 
distortion, photometry, and cosmic ray cal- 
ibrations). For each image, we create a 
map of pixels that deviate by more than 
3.5 u in a 20x20 pixel box and we replace 
them with the median value in the box. We 
multiply each image by the photflam of 
its header to convert the pixel counts to 
erg s _1 cm -2 A -1 arcsec -2 . For the Fo- 
malhaut PSF registration, we first start 
with the 2006 sequence at F606W that is 
recorded with the 3.0" focal plane mask. 
We align every image maximizing its corre- 
lation with the first image of the sequence 
(Table 1). The correlation is maximized in 
the annulus with inner and outer radii of 
140 and 200 pixels where the central ver- 
tical band of width 240 pixel is removed 
(saturated detector) and where the coro- 
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Fig. 1. — ACS images of the dust belt and the object b (circle) around Fomalhaut at F606W 
in 2004 (top) and 2006 (bottom). Two arrows point to background sources. The length of 
the segments giving the East and North orientations is 2.5" . The intensity scale is linear 
and it is the same for the two images. 
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Fig. 2.— Same as Fig.Ulfor F435W (top) and F814W (bottom) ACS images taken in 2006. 
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Fig. 3. — Same as Fig. {T\ for the STIS image taken in 2010 with the CLEAR filter. 



nagraphic focal plane bar is masked. We 
call this optimization area A. Once the im- 
ages are aligned to within 0.1 pixel, the ab- 
solute center of the star PSF is found by 
median-combining the aligned images and 
by registering the resulting image to the 
image center by maximizing in A the cross- 
correlation of the diffraction spikes with 
themselves in a 180° rotation of the im- 
age about its center. This procedure de- 
fines the absolute center within 0.5 pixel 
and we call R606 the registered median- 
image. We then register all images of 
F606W sequences maximizing their cross- 
correlation with R606 in A. For the F435W 
and F814W sequences, we scale R606 to 
the corresponding wavelengths (814nm for 
F814W and 480nnfl for F435W) and call 



R435 and R814 the resulting images. We 
then register all F435W and F814W im- 
ages maximizing the cross-correlation in A 
with R435 and R814 respectively. For ev- 
ery sequence listed in Tab.Q], we then sub- 
tract the stellar speckles. As the field 
of view rotates only by a few degrees, 
there is almost no difference between ap- 
plying a locally op timized combination of 
images algor i thm (iLafreniere et al.l 120071 ; 



Marois et al.l l2010bl ) or a basic angular 



1 Better match of the diffraction spikes than for 
435nm. 



differentia l imaging data redu ction as de- 
scribed in lMarois et al.l ( 2006 ). We choose 
the second procedure which is less time 
consuming. Considering one of the se- 
quences, we subtract from each image a 
reference PSF that is the median of all im- 
ages that were recorded at a different roll 
angle. We then rotate the images to put 
North up and median combine them. For 
the 2006 data, we work out the weighted 
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mean of the reduced images taken with 
the 1.8" and 3.0" masks in the same fil- 
ter. As the fields-of-view do not exactly 
overlap, the contrast is not the same in 
all parts of the images ( Figs. [1] and [2]). A 
point source (arrow) is detected in all im- 
ages South West from Fomalhaut A. An 
extended object (red arrow, South East) is 
also detected in 2004 and in 2006 (F814W). 
The motion of these two sources between 
2004 and 2006 are consistent with them be- 
ing background objects. Fomalhaut b (in- 
side the circles) is detected at F435W, 
F606W, and F814W with a signal-to- noise 
ratio of ~ 3, 5-6, and 3 respectively and it 
does not have the same motion as the back- 
ground sources (Fig.HJ). We thus confirm 
that it is a real object orbiting Fomalhaut. 



Fomalhaut 



10.2 



10.0 



Expected bckgd 
Measured at t 
■ Measured 




-9.0 -9.5 -10.0 

RA (arcsec from star) 



-10.5 



Fig. 4. — Fomalhaut b measured positions 
in 2004, 2006 and 2010 images (red crosses) 
and expected positions for a background 
source (green crosses). 



3.2. STIS 

The sx2 images that are provided by 
the STIS pipeline (geometric distortions, 
photometry, and cosmic ray calibrations) 
are used for our analysis. The flux density 
is converted to erg s _1 cm" 2 A" 1 arcsec -2 
by multiplying each image by the phot- 
flam of its header and dividing it by the 
exposure time. The spider spikes are well 
detected in these images and we register 
the first image of the sequence maximiz- 
ing the cross-correlation of the spikes with 
themselves in a 180° rotation of the im- 
age about its center. The maximization 
was done around the spikes (±2 pixels) be- 
tween 100 and 230 pixels from the star. 
We then register the other images maxi- 
mizing the cross-correlation with the first 
image in the 200 pixel-radius disk where 
the 160 pixel central vertical band and 
the 30 pixel central horizontal band are 
masked. As the roll angles are well spread 
in the 0-157° interval, we apply a locally 
optimi zed combination of images algo- 



rithm ( Lafreniere et al.ll2007l ; iMarois et al. 



2010bl ) to suppress the stellar diffraction 
pattern. Using a LOCI algorithm, a PSF 
reference image is built for each image of 
the sequence and it is subtracted. After 
subtraction, the images are rotated to put 
North up and they are median-combined. 
The final image with the detection of Fo- 
malhaut b is shown in Fig. El 

4. Data analysis 
4.1. Belt geometry 



have been discussed previouslv ( 


Kalas et al. 


2005; 


Acke et al. 


2012; 


Bolev et al. 


2012) 



and it is beyond the scope of the present 
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paper to refine them. However, we find an 
eccentric belt that reproduces the images, 
with an eccentricity e=0.10-0.11, a radius 
between 136-148 AU, a right ascension of 
ascending node 156.5-157°, an argument 
of periapsis 35°, and an inclination of 67°. 
All the parameters are in good agreement 
with the published values of K08 (31±6° 
for the argum ent of peri a psis u nlike the 
1±6° found bv lAcke et all faoiti )). We did 
not use a mathematical fit to optimize val- 
ues of parameters and our best visual fit is 
only used to estimate the belt geometry in 
our images. 

4.2. Astrometry 



the difference Si between K08 positions and 
ours at epoch i as 



We use the Tiny Tim (IKrist et al.ll201ll ) 
tool that generates HST template PSFs 
to build a model of a point source in 
our images at the position of Fomal- 
haut b to accurately estimate its astrom- 
etry and photometry. The temperature 
of a thermally emitting companion that 
would reproduce the F814W flux is roughly 
400 K (j.Tanson et al.l I2012L hereafter J12). 
Thus, we choose a 400 K source in the Tiny 
Tim tool and we simulate the images prior 
to the speckle suppression to account for 
the ADI/LOCI effects, the rotation we ap- 
ply to put north up, and the weights of the 
weighted means for the 2006 data. We then 
adjust the position and flux of the tem- 
plate to subtract from the image to min- 
imize the residual noise in a 0.7" (1.2")- 
diameter disk centered on Fomalhaut b 
for the ACS (SITS) data. The positions 
we derive from the fit are given in Tab. [21 
K08 measured in their images that Foma- 
lhaut b is at [ra, dec] = [-8.62", 9.20"] 
and [-8.60", 9.38"] from Fomalhaut A in 
2004 and 2006 respectively. We estimate 



' 5jra 2 



5,- dec 



a 



dec,i 



where <5jra and <5jdec are the difference 
between K08 measurements and ours of 
the offset along the West-East direction 
and the South-North direction respec- 
tively. <r rai i and o"dec,i are our error bars (K08 
give no error bars). We find that our po- 
sitions are within 1.5 a of K08 positions 
at the two epochs 2004 and 2006 (i.e. 
5i < 1.5). The difference with K08 could 
result from a different registering technique 
or from differences in the ACS pipeline 
that have been upgraded since 2008. 

As we have only three epochs close in 
time and large error bars for the 2010 data, 
we cannot strongly constrain the orbital 
parameters. We then consider only two ke- 
plerian orbits - one that crosses the dust 
ring and a second that does not cross the 
ring - and compare the expected and mea- 
sured positions. 

The first orbit is a 0.19 eccentric or- 
bit with a 118 AU semi-major axis, a 156° 
right ascension of the ascending node, a 70° 
inclination, and a 2° argument of periap- 
sis. The orbit does not cross the dust 
ring and is represented in dashed lines in 
Fig. [5] where the dust belt is bound by 
dashed-dotted lines. We use Eq.[TJ replac- 
ing K08 positions with the expected posi- 
tions of Fomalhaut b on the keplerian orbit 
to estimate the differences Si between the 
expected positions and our measurements 
at each epoch % (2004, 2006, and 2010). 
Then, we estimate the total difference as 
S = J2i We find that the expected po- 
sitions are 1.1 a from the measured posi- 
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Date 




Filter 


Offset from A to b 
RA (") DEC (") 


Separation 

(") 


PA 

(deg) 


2004 Oct. 


25-26 


F606W 


-8.59±0.02 


9.19±0.02 


12.58±0.03 


316.9±0.1 


2006 Jul. 


14-15-16 


F435W 


-8.61±0.03 


9.36±0.03 


12.72±0.04 


317.4±0.2 


2006 Jul. 


19-20 


F606W 


-8.64±0.02 


9.36±0.02 


12.73±0.03 


317.3±0.1 


2006 Jul. 


18-19 


F814W 


-8.64±0.03 


9.36±0.03 


12.73±0.04 


317.3±0.2 


2010 Jun. 14-Sep. 13 


CLEAR 


-8.81±0.07 


9.79±0.07 


13.17±0.10 


318.0±0.4 



Table 2: Fomalhaut b SLstrometry with respect to Fomalhaut A. 



tions (5 = 1.1). The second keplerian or- 
bit (full lines) we consider has an eccentric- 
ity 0.25, a semi-major axis 119 AU, a right 
ascension of the ascending node 156°, an 
inclination 70°, and an argument of peri- 
apsis 10°. The difference between expected 
and measured positions is 0.9<r. If Foma- 
lhaut b follows this orbit, it will enter the 
dust belt in ~ 100 years. Given these two 
examples, new data are obviously required 
to conclude whether Fomalhaut b will or 
will not cross the belt. 

4.3. Fomalhaut b as a point source 

We consider Fomalhaut b as a point 
source in this section. We estimate its pho- 
tometry and compare our results with K08 
fluxes (§ !4.3.ip . We then examine mod- 
els discussed by K08 and J12 (§133 
and § ESD - 

4-3.1. Photometry 

For each filter and epoch, we derive the 
photometry by integrating the flux den- 
sity of the PSF template that best fits the 
data (§ 14.21) in a 1.25" radius disk for the 
ACS data and 1" radius for the STIS data. 
The fractions of the PSF integrated energy 
inside these apertures are 0.960, 0.961, 
0.918, and 0.996 for the F435W, F606W, 



F814W, and CLEAR/STIS filters respec- 
tively. To convert the estimated flux den- 
sities Fx in erg s _1 cm -2 A -1 to flux densi- 
ties F v in erg s -1 cm -2 Hz -1 (i.e. 10 23 Jy), 
we use the photplam keyword recorded 
by the ACS and STIS pipelines in the fits 
headers: 

F v = Fx photplam 2 10" 18 - 4768 ( 2) 

The resulting flux densities (/iJy) are given 
in Tab. [21 The error bars o v in percent- 



Date 


Filter 


Flux density (/iJy) 








Kalas 


2004 


F606W 


0.69±0.11 


0.61±0.05 


2006 


F435W 


0.41±0.15 


<0.87 (5<r) 


2006 


F606W 


0.47±0.09 


0.29±0.03 


2006 


F814W 


0.38±0.12 


0.37±0.04 


2010 


CLEAR 


0.88±0.32 





Table 3: Photometry if Fomalhaut b is a 
point source. 



age are the inverse of the signal-to-noise 
ratios. In these ratios, the signal is the 
integrated flux density inside a 0.2" radius 
disk centered on Fomalhaut b and the noise 
is the square root of the total variance of 
the residual noise after subtraction of the 
best PSF in the same area. K08 express 
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Fig. 5. — Trajectories that best fits the Fomalhaut b positions with eccentricity 0.19 (dashed 
lines) or 0.25 (full line). 
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their Fomalhaut b photometry and upper 
limits in Vega magnitudes. We convert 
their measurements to /iJy (last column 
in Tab.E]) using the ACS handbook. 

Most of the flux densities (Tab. [3]) 
are consistent with K08 values except 
our F606W/2006 point which is ~ la 
brighter [a is the quadratic sum of K08 
error bars and ours). Moreover, our er- 
ror bars are larger than K08's ones. Thus, 
even if we still detect a variability in the 
F606W filter between 2004 and 2006, it is 
not as significant (1.5 a a) as it is in K08 (5- 
6 ok) - where o~g and o~k are our error 
bars and K08 ones respectively. We also 
find that the flux density measured in the 
CLEAR filter (its bandpass roughly cor- 
responds to F435W+F606W+F814W) is 
consistent with the three ACS flux den- 
sities given the large error bar. Finally, 
we (marginally) detect Fomalhaut b at 
F435W unlike K08 who have an upper 
limit. 

We plot the photometry of our de- 
tections (crosses) in Fig. [6] along with 
upper limits from the literature (K08, 



Unresolved Fomalhaut b 



Marengo et al.l 12009k and J 12) at vari- 



ous wavelengths. J12 find that to comply 
with their 4.5 /mi upper limit, the plane- 
tary mass upper limit is lMj at 400 Myr. 
Thus, we compare the measurements with 
a model of a cloud-free atmosphere for 
a 1 Mj planet at 400 Myr with th e so- 
lar metallicity (ISiegelfeBurrowsll2012l ) (full 
line in Fig. [6]). It is clear that a planet-only 
model for Fomalhaut b is not consistent 
with the observations. K08 proposed two 
other models: a cloud of dust (§ I4.3.2|) or a 
disk of dust around a Jovian-planet (§ 14.3.31) . 
We revisit these two models in light of our 
updated photometry. 



1000.0 



100.0 



°> 10.0 




Wavelength (/im) 

Fig. 6. — Fomalhaut b flux density ((iJy) 
for various wavelengths (fim^ 1 ) in the case 
the object is unresolved. Crosses corre- 
spond to our detections (four black for 
ACS and one light blue for STIS). Ar- 
rows are upper limits from the literature: 
green, red, an d blue for K08, J 12, and 
Marenso et "ail feOQSl ) respectively. The 
solid line represent a cloud-free atmosphere 
model for a 1 Mj planet at 400 Myr. 

4.3.2. Cloud of dust 

We consider the model introduced in K08 
with a 0.53 AU diameter cloud composed 
of dust grains with a differential size dis- 
tribution dn/da oc (a/a ) -3 ' 5 where the 
radius a goes from a min to 1000 /im. Using 
Mie theory, K08 calculate the apparent 
magnitudes of such a cloud composed of 
water ice (density=l, mi cc ) or refractory 
carbonaceous material (density=2.2, ttilg) 
with a m ; n = 0.01 /im (hereafter m 001 ) 
or 8 /tm (hereafter m 8 ). The total mass 
in grains is adjusted such that the inte- 
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grated light in F814W from the model 
matches K08's observations (K08's and 
our photometry in F814W are in agree- 
ment). We convert the Vega magnitudes 
provided in K08's Tab. S3 to flux densities 
in /iJy (Tab.HJ). The last line gives the 
error e between the expected flux densi- 
ties F BjU and the observed densities F v : 



e 




(F v - F e ,„) 2 


( 


Filter 




mO.Ol 
LG 




m LG 


F435W 
F606W 
F814W 


0.71 
0.55 
0.37 


0.58 
0.50 
0.37 


0.63 
0.46 
0.37 


0.46 
0.45 
0.37 


e 


2.0 


1.3 


1.7 


1.0 



Table 4: Expected flux densities F e ^ (fj,Jy) 
derived from the cloud models (m^ 1 , 
m LG> m fce; an d m LG) proposed in K08. 
The last line gives the difference be- 
tween F E)i , and our measured photome- 
try F v (see text for details). 



K08 reject the possibility that Fomal- 
haut b can be explained by one of these 
cloud models because 1/ they do not detect 
the object at F435W (they do not reject 
?til G for this reason), 2/ the red color they 
observe does not match the model, and 3/ 
they cannot explain the F606W variabil- 
ity. All these reasons do not apply to our 
new photometry because 1/ we detect Fo- 
malhaut b at F435W, 2/ the expected flux 
densities match the observed flux densi- 
ties within 1.7cr for three of the four mod- 
els (e < 1.7), and 3/ the F606W variability 
is not significant in our images. K08 also 
explain that such a cloud could result from 



a collision of two planetesimals and that 
the probability of such an event is lower at 
the Fomalhaut b position than closer to the 
star or closer to the belt. However, as sug- 
gested by J12, the probability of a collision 
is not the probability of its detection be- 
cause the speckle noise and the high bright- 
ness of the ring may prevent detections of 
such clouds close to the star and the belt 
respectively. Finally, K08 argue that such 
dust clouds would be sheared due to differ- 
ential gravitational forces and rapidly spa- 
tially resolved by HST. However, assum- 
ing a cloud with diameter 0.5 AU only sub- 
ject to gravitational forces from the star, 
we find that its image would be larger 
than 2 pixels (~fwhm) and 4 pixels after 
~100 years and ~200 years respectively. It 
would take ~500 years to shear the cloud 
of dust so that it could be spatially re- 
solved in the HST images with no doubt. 
Thus, we find no strong arguments to re- 
ject K08 models of a dust cloud with ra- 
dius ~ 0.5 AU, composed of water ice or 
refractory carbonaceous small grains, and 
younger than ~ 500 years. 

4-3.3. Material surrounding a Jupiter- 
like planet 

A second scenario proposed by K08 is 
an unseen Jovian planet surrounded by a 
disk of dust with a radius of 16-35 planet 
radius. As the K08 photometry is close to 
ours and K08 only work out rough numbers 
(they could not constrain all the param- 
eters with only two photometric points), 
the 16-35 planet radius disk surrounding 
an undetected Jupiter-like planet is con- 
sistent with our photometry. J12 reject 
this model because 1/ it does not explain 
the F606W variability and 2/ the belt ge- 
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ometry would be strongly affected con- 
sidering a ring-crossing o rbit for Fomal- 
haut b ( Kalas et al. 2010l ). As we do not 
find a significant F606W variability and 
our new astro metry cannot reject an orbit 
that does not cross the ring, we cannot rule 
out this model using the K08 arguments. 
J12 also consider that if the spin of the 
star is aligned with the plane of the disk, 
the north-west side of the disk is closer to 
Earth than the south-east side. In that 
case, Fomalhaut b is between its star and 
Earth in the radial direction and J12 claim 
that it would be difficult to explain how 
an optically thick disk can reflect so much 
light towards Earth. It is true if we ob- 
serve the non-illuminated side of the disk 
but we can imagine an inclined disk such 
that we observe the illuminated part of the 
disk even if Fomalhaut b stands between 
Fomalhaut and Earth. 

4.4. Is Fomalhaut b resolved ? 

4-4- 1- Extended source vs PSF 

Given that a possible model for Foma- 
lhaut b involves a cloud of dust, it would 
be possible that Fomalhaut b has slowly 
expanded in time and we test here the pos- 
sibility that the Fomalhaut b HST images 
are slightly spatially resolved. To do so, we 
consider the same PSF templates as in § 14.21 
and we model the object intensity distri- 
bution as a uniform disk with radius R. 
We convolve the object model by the PSF 
template and obtain the object image tem- 
plate T . We then calculate a signal-to- 
noise ratio SNR. The signal S is the in- 
tegrated energy of the observed Fomal- 
haut b inside a disk of radius R+fwhm/2. 
The noise N is the standard deviation of 
the residual noise inside a 12 pixel radius 



disk around Fomalhaut b after subtraction 
of T Q . As the radius of 12 pixels is larger 
than R+fwhm/2, N is large if the mod- 
eled source is too large or too small to fit 
the data. As S and N are not calculated 
over the same surface, we normalize N by 
the square root of the ratio of the number 
of pixels in the 12 pixel radius disk to the 
number of pixels in the R+fwhm/2 radius 
disk. Finally, SNR is the ratio S/N. 

For each radius R, we search for the 
flux to apply to T Q to maximize SNR. The 
center column in Figs. [7] and [8] give the 
evolution of SNR with the radius R in the 
2004/F606W, 2006/F435W, 2006/F606W, 
2006/F814W, and 2010/CLEAR images 
from top to bottom (y-scales and x-scales 
are not the same for all plots). The images 
of Fomalhaut b and of the residual noise 
after the best subtraction are overplotted. 
Note that even for the STIS data (last 
row), we express the radius in ACS pixels, 
i.e. 25 mas. 

In the 2006/F814W filter (center of the 
fourth row), the SNR criterion is maximum 
for an extended source with R=5 pixels 
(before convolution). Fig.0 presents the 
residual intensities after subtracting a 
PSF (center) or an extended source (right) 
from the Fomalhaut b image (left) at 
F814W. We expect to find the same size in 
all images but the source does not seem to 
be as large in the other filters (Figs. [7] and 
E}. This may be due to biases in our SNR 
criterion because of the low signal-to-noise 
ratio of the detections. To estimate these 
biases, for each filter/epoch, we add fake 
sources to six speckle noise dominated ar- 
eas in the real images (different colors in 
Figs. [7J and [8]) with the signal-to-noise level 
of Fomalhaut b and we apply the same 
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Other P5F F606W 



Other Ext F606W 




Extended source radius 



Fig. 7. — SNR criterion as a function of the source radius in ACS pixels (25mas) for a 
synthetic point source (left column), the Fomalhaut b image (center column), and a synthetic 
extended source (right column). Top to bottom, the rows correspond to F606W in 2004, 
and F435W and F606W in 2006 respectively. See text for details. 



procedure to plot SNR function of R. 
We consider two different fake sources. 
We first assume an unresolved source (left 
column) and we simulate its image us- 



ing 25 TinyTim PSFs at different posi- 
tions on the detector (one curve per PSF). 
We do so to ensure that the PSF vari- 
ations cannot mimic an extended source 



14 



Other PSF FS14W 



Other Ext F814W 




Extended source radius in pixels 
Other P5F STIS 



Extended source radius in pixels 
2010 STIS 



FWHM - 2.0 




FWHM - 2.0 




8 10 



Extended source radius in pixels 
Other ext STIS 



FWHM - 2.0 " 


















l Fake*V 







Fig. 8.— Same as Fig.\7\for F814W in 2006 (top) and CLEAR in 2011 (bottom). 
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Fig. 9. — Fomalhaut b at F814W (left) after subtracting a PSF (center) or an extended 
source (right). 



for our SNR criterion. The second fake 
source has an intensity distribution uni- 
form in a disk with R=3 pixels before con- 



volving with the PSF (i.e. 1.16 AU diam- 
eter, right column). We use the same 25 
TinyTim PSFs to convolve the object and 
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obtain 25 images. For the STIS data, only 
one fake PSF is used because the TinyTim 
tool does not account for the PSF varia- 
tions with the positioning on the detector. 
For the 2006/F606W (third row), the ob- 
servations (center) can be reproduced by 
both the point source (left) or the extended 
source (right). The curves of the other 
observations seem closer to the extended 
source case. 

In summary, we find that the Foma- 
lhaut b HST images are better repro- 
duced by a slightly extended source and 
we choose in the rest of this section a basic 
model of a uniform disk with a diameter 
of 1.16 AU before convolution by the PSF. 
However, we are aware that the signal-to- 
noise ratios of the detections are low and 
that new observations are required to fully 
confirm this result. 

4-4-2. Unlikely an instrumental effect 

We investigate what instrumental effect 
or data processing could explain such an 
extended source. 

The ACS/HRC PSF is contaminated 
by a halo for red sources, especially at 
F814W. The halo which adds to the "nor- 
mal" PSF has a diameter (42-2.36 A) pixels 
and contains a total fractional intensity 
2 (A — 0.45) 3 for the wavelength A in mi- 
crons. A 10 pixel diameter halo requires 
a dominant flux at A ~ 1 1 urn. from this 
expression, which does not make sense be- 
cause it is well outside the sensitive band- 
pass of the detector. Moreover, even if the 
signal-to-noise ratio is low, we do not ob- 
serve in the F814W image a PSF plus a 
halo but only an extended image. 

A second explanation for such an im- 
age could be a misregistering of the raw 



images. In that case, after the rotations 
that put north up in the ADI process, 
all the Fomalhaut b images would not fall 
at the exact same position, resulting in a 
blurred image. If this happens, any source 
in the field of view would be affected the 
same way. As the background source image 
(south-west) is very similar to the Tiny- 
Tim PSF for the three filters (Fig. fTOl) . we 
can rule out this possibility. 




Fig. 10.— 2006 F814W image (left) and 
Tiny Tim model (right) of the south-west 
background source. 

The last instrumental effect that we 
foresee is a differential geometric distortion 
of > 1 pixel at the Fomalhaut b position 
between the images of a same sequence. 
The ACS pipeline corrects for the distor- 
tions with an accuracy 0.01 pixels (ACS 
handbook[§). Thus, it would require dif- 
ferential distortions 100 times larger than 
the pipeline accuracy at the Fomalhaut b 
position but almost no distortions at the 
background source position (Fig. fTOj) which 
is roughly at the same angular separation 
from the star. This scenario seems very 
unlikely. 

Finally, we find no instrumental effects 
that could explain the possible spatial res- 
olution of Fomalhaut b in our images. 



2 http: / / www.stsci.edu/hst / acs / documents /handbooks /cycle20 /c!0 Jr 
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4-4- 3. Photometry 

For each filter/epoch, we consider the 
template T for a 1.16 AU diameter object. 
We adjust its flux to minimize the resid- 
ual noise in a 0.7" radius disk when we 
subtract it from the observations. For the 
STIS data, the disk radius for the mini- 
mization is 1.2". We follow the steps de- 
scribed in § 14.3.11 to convert the flux densi- 
ties to Jy and estimate the error bars. The 
results are given in Tab.|5j As expected, 

Date Filter Flux density (/iJy) 



Unresolved Fomalhaut b 



2004 F606W 

2006 F435W 

2006 F606W 

2006 F814W 

2010 CLEAR 



1.09±0.13 
0.71±0.14 
0.67±0.11 
0.62±0.12 
1.26±0.34 



Table 5: Photometry if Fomalhaut b 
is 1.16 AU large. 

the fluxes are larger than in the case of 
a point source. Moreover, the flux varia- 
tion at F606W is larger than in the point 
source case but it is still less than 2.5 aa, 
thus not yet significant. An unfortunately 
situated speckle at less than 3 pixels from 
Fomalhaut b could explain this variation. 
Finally, the flux density measured in the 
large band of STIS is consistent with the 
average flux density measured in the ACS 
filters within 1.4 o~q. In the case we resolve 
Fomalhaut b, we plot the photometry of 
our detections (crosses) in Fig.fTTl 

The models presented in sections 14.3.21 
and 14.3.31 of a planet surrounded by a dust 
disk and a cloud of dust are both intro- 
duced in K08 in the case Fomalhaut b is 
not resolved. In the case of a spatially re- 
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Fig. 11. — Same as Fig.\B[ in the case the 
object is resolved. 

solved Fomalhaut b, we propose one basic 
model that assume that Fomalhaut b is the 
result of the collision of two Kuiper belt ob- 
jects (§ !4.4.4p and we adapt a model of cir- 
cumplanet ar y swarm of satelli t es (§14.4.51) 



proposed by iKennedyfcWyattl (1201 IT ) 



4-4-4- Collision of Kuiper belt objects 

In this section, we propose a basic 
model to roughly estimate the size and 
the amount of light that is scattered by a 
cloud of dust produced by the collision of 
two Kuiper belt objects (KBO). The objec- 
tive is not to derive the exact radius, mass 
and velocity of the KBOs that could create 
Fomalhaut b but to show that the collision 
of two KBOs is not completely inconsistent 
with the observations. First, we estimate 
the total grain mass that can explain the 
fluxes received from Fomalhaut b. Then, 
we show that the amount of dust can be 
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the result of a collision of two 50 km radius 
colliders. We evaluate the rate of collisions 
of two such KBOs around Fomalhaut. Fi- 
nally, we estimate when the collision may 
have occurred to reproduce the size of the 
Fomalhaut b images. 

If the particles of dust are spheres with 
radius a and if the cross section of the par- 
ticles equals their geometric albedo, the 
mass Met of a cloud of dust t hat lies at 



a dis tance D from its star is (jJura et al. 
1995h 

M rf >^^ 2 a| (4) 

where p is the mass density of the dust 
grains, and L sc and L* are the luminos- 
ity of the light scattered by the cloud and 
the luminosity of the star respectively. In- 
stead of estimating the ratio of the lu- 
minosities, we work with the fluxes F sc 
and received at the telescope. We es- 
timate F sc from the measured fluxes (F) 



in the ACS filters (i= 
F814W, Tab. ED 



:F435W, F606W, and 



F ^ 



(5) 



with Az/j the bandwidths of the filters. Our 
estimation of F sc only includes the scat- 
tered energy in the F435W, F606W and 
F814W bandpasses. F* has to be calcu- 
lated for the same bandpass. Assuming a 
Planck law, F* in the bandpass [\ m in = 



435 - 50 nm, A r 



825 + 115 nm is 



rp J Umin 1 v ' fc\ 

J u I (exp u — 1) dti 



where u min , max = hc/(kT \ max>min ) with 
the Planck constant h, the speed of light 
in vacuum c, the Boltzmann constant k, 



the stellar effective t emperature T (8751 K, 
Di Folco et al 1 l2004h . and the stellar flux F* tot 
received at the tele scope (8.914e-6erg.cm~ 2 , 



Kalas et all 120081 ). For D ~ 120 AU 
and dust grains with radius a =10 /im 
and p = 2g.cm -3 , we find from Eqs.HJ 
and |6] that the total grain mass needed to 
reproduce the photometry of Fomalhaut b 
is M d ~ 4.10 19 k. 



Jewittl ( 120121 ) estimates the mass m c of 



particles that are ejected after a collision of 
two KBOs with a mass M kbo , a radius r, a 
density p, and a relative velocity U 



Mi 



A 



kbo 



SnGp 



n -1.5 



u 



1.5 



(7) 



6.6710- n m 3 



where A equals 0.01 and G 
kg" 1 s~ 2 is the gravitational constant. 



Jewittl (120121 ) assumes the particles have 
radii in the range 0.1 pm < a < 0.1 m with 
a power law distribu tion in radii with in - 
dex ~ 3.5 (see also iKadono et al.l |2010| ). 



For typical KBOs in the ring, U is the or- 
bital velocity times h r /(2D) with h r the 
full vertical height of the ring at radius D. 
With h r ~ 3.5 AU (jKalas et al. l 120051 ) at 
D ~ 120 AU, U is close to 60 m s -1 . As- 
suming KBOs with radius r = 50 km, the 
total debris mass m e after the collision is 
roughly 1% of the mass M kbo of one of the 
two colliders with a density p = 2 gem -3 . 
Given the approximations in the models, 
the expected mass of dust (l%Mkbo) that 
is ejected after a collision of two 50 km 
radius KBOs is consistent with the mass 
estimated from the photometry of Fomal- 
haut b (4 % M kbo for a 50 km radius KBO). 

If the escape velocity r y/8 ix G p/3 is the 
average velocity of the particles, the diam- 
eter s of the cloud is 2r t^/8nGp/3 at 
the date t after the collision. To reproduce 
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the observed size s = 1.16 AU (§ 14.4. ip . the 
collision should have occurred ~ 50 years 
ago. We can also estimate from this ex- 
pression that after ~ 200 years, the source 
would have a diameter ~ 4AU and would 
be ~ 10 times fainter than the current de- 
tections assuming the same amount of re- 
flecting dust. It would not be detected in 
our images. As the size of the possible 
extended source is very approximative, we 
keep in mind that these numbers are coarse 
estimations. 

Finally, we evaluate the rate of a colli- 
sion of two r = 50 km KBOs at the Fomal- 
haut b position. We first find the collision 
time which reads 



1 



Airnr 2 U 



with n the number of KBOs with ra- 
dius r per unit volume. To estimate n 
we need the mass Mdi S k of the debris 
disk around Fomalhaut. We assume 
Mdisk=40 Msarth because 1/ estimates of 
the mass of the Sun's early Kuipe r belt 
is 40M Ea rth (ISchlichtingfcSaril bOQfil ). and 
2/ estimates of the mass of the Vega de- 
bris belt i s 10 Mparth in objec ts with radii 
< 100 km (jMuller et al.ll2010l ) whereas the 
Vega IR luminosity is 4 times fainter than 
that of Fomalhaut. Simplifying by setting 
the radii of all the KBOs to 50 km does not 
qualitatively alter the collision frequency 
estimated below. Under these conditions 
and considering the belt surrounding Fo- 
malhaut A has a yolume 2 n D AD h r with 
AD ~ 0.13 D dKalas et al l l2005h . the 
number of KBOs is 

n 2 7T D AD h r 



M-i; 



disk 



Mi 



2 x W 



kbo 



(9) 



Using U = Txh^MjM Q {lK\]/DY, and 



Eqs.[8] and[9l we can write 



0.13 f DY Mi 



2tt 



'kbo 



M 



disk 



D 
1AU 



M m 



M* 
(10) 

Finally, the rate k of collisions of two 
KBOs with radii r is the ratio of the num- 
ber of KBOs (Eq.ED to t col 



H 



2tt ( r y (M Ais] 



0.13 \DJ V M 




1AUV M* 



(kbo/ v v D J M 
(11) 

Eq.EUwith M disk = 40M Earth , and M* = 
2 M Q indicates that ~ 1 collision of two 
r = 50 km KBOs occurs every century 
in the ring around Fomalhaut A. Even if 
the rate at the Fomalhaut b position may 
be lower, given the approximations of our 
model the rate is high enough to make 
such a ~ 50 year-old event plausible. At 
the same time, it is low enough to explain 
that we detect only one event around Fo- 
malhaut as each event would be detectable 
during ~ 200 years in our images. 

In summary, we conclude that it is plau- 
sible that Fomalhaut b is a cloud of dust 
that was produced ~ 50 years ago by the 
collision of two KBOs with radii ~ 50 km. 

4-4-5. Circumplanetary satellite swarm 



KennedvfeWvattl (120111 . Kll) propose a 



model of circumplanetary satellite swarms 
that they apply to Fomalhaut b. They find 
that the planet mass can be ~2-100 Msarth 
surrounded by a swarm that lies at 0.1-0.4 
Hill radii. The swarm mass would be of 
the order of a few lunar masses. But these 
numbers are derived from K08 photometry 
of an unresolved source. 

Here, we use the same model under 
the same assumptions (body size dis- 
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tribution and maximum/minimum body 
sizes, dust density, etc) but we consider a 
swarm of satellites with diameter 1.16 AU, 
our photometry (Tab. EJ, and a star with 
age 4 40 Myr instead of 200 Myr (jMamajek 
20121 ). We do not describe the model as it 
is done in Kll. We only use the meaning- 
ful equations to constrain the planet mass 
and the swarm mass and size following the 
steps in §3.3.1 and §3.3.2 of Kll. First, 
we derive the total cross-sectional area of 
dust cr tot from our photometry: <r tot = 
7.22 x 10~ 4 au 2 , assuming a geometric 
albedo 0.08, a phase function 0.32 (Lam- 
bert sphere at maximum extension from 
its host star), the star effective t empera- 
ture 8,751 K flPi Folco et al.ll2004h . and a 
stellar luminosity 6.34.10 27 W (K08). 

As we assume that we resolve Fomal- 
haut b, we can write 2 77 Rmn = s, where 77 
is the semimajor axis of the satellites of 
the swarm relative to the Hill radius Rmn 
at the Fomalhaut b separation (118 AU) 
and s is the swarm diameter. As s = 
1.16 AU is an approximative size (§ 14.4. 1 j) . 
we consider 1.16 AU< s <2.32AU. Us- 
ing the Hill radius expression and 2 solar 
masses for Fomalhaut (Kll), we derive two 
constraints: 0.62/Mp 1 / 3 < r] < 1.24/Mp 1 / 3 , 
where M p \ is the planet mass. 

If the cloud is optically thick and the 
dust is uniformly projected on the sky 
across n (77 Rmn) 2 , the geometrical opti- 
cal depth r = o- t ot/(^(v -Rmn) 2 ) is smaller 
than 1 (Kll), which sets r] > 0.017/M p Y 3 . 

Then, considering a collision-limited 
satellite swarm around Fomalhaut b (i.e. 
swarm has just started to suffer colli- 
sions) that reproduces the observed cr to t, 
it imposes a minimum limit for the satel- 
lite semimajor axis rj > 0.33/Mp| 12 for 



a 440 Myr system (see Kll for details). 

Kll also study the collision velocities 
that are required to destroy a large ob- 
ject at the Fomalhaut b position. Us- 
ing their expressions in the case of a re- 
solved object we set a constraint that reads 
rj > 0.71/M°- 46 (Kll). 

Moreover, Kll assumed that satel- 
lite orbits with 77 > 0.5 are not sta- 
ble and did not consider them. Finally, 
we account for the 1 Mj upper limit 



that iJanson et al.l (120121 ) put from the 
non detection at 4.5/im for a 400 Myr sys- 
tem (clo se enough to 44 ± 40 Myr pro- 
posed by lMamajekll2012l ). We plot all the 
constraints in Fig. [12] that gives the semi- 
major axis 77 of the satellites against the 
planetary mass M p \. The parameters for 
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Fig. 12. — Semimajor axis rj of satellites of 
the swarm versus planetary mass M p \ dia- 
gram that shows the different constraints 
derived in the text. The parameters that 
could explain the Fomalhaut b images are 
inside the dashed region. 
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which Kll's model can reproduce the pho- 
tometry of a 1.16 AU source are within the 
dashed area. The minimum and maximum 
planetary masses are ~2 Msarth and 1 Mj 
and the swarm has a total mass 3-12 M Mo on 
and lies at 0.17-0.5 Hill radii around the 
planet. 

In the case of an unresolved object (Fig. 7 
in Kll), Kll find that the mass of the 
planet (< 100 M Earth ) is not sufficient for 
it to have a significant gaseous envelope 
and enable mechanisms that could explain 
the migration of Fomalhaut b that presum- 
ably originates somewhere closer to the 
star. Kll also argue that a single planet 
with mass < 100 MEarth - which is simi- 
lar or less than the mass of the main de- 



bris ring (1-300 M F , aTth , IWvattfeDentll2002 



Chiang et al.l 120091 ) - is unlikely responsi- 
ble for shaping the dust belt. In our case of 
a source with diameter 1.16 AU, the range 
of the planetary mass goes up to 1 Mj and 
a Jupiter-like planet can have a significant 
gaseous envelope and shape the dust belt. 

5. Conclusions 

Our independent analysis of the ACS 
and STIS data taken in 2004, 2006, and 
2010 confirms that Fomalhaut b is real and 
is not a speckle artifact as we clearly detect 
the object at the three epochs at several fil- 
ters (Fi gs-tHEl and[3]). In t his way, we con- 
firm the iKalas etloTfeoOSl . K08) detection. 
However, we find differences in our analy- 
sis concerning astrometry and photometry 
of Fomalhaut b. 



Unlike IKalas et al.l ( 120101 ) , we cannot af- 
firm that the object follows a trajectory 
that crosses the belt of dust because our 
astrometry is consistent within 1.1 o with 
crossing and non-crossing orbits (§ 14. 2\ . 



We detect Fomalhaut b in the short 
wavelength filter F435W whereas K08 find 
an upper limit. In the case of an unresolved 
source, our photometry is consistent with 
K08 at F606W/2004 and F814W/2006 but 
differs at F606W in the 2006 data (§SXTJ. 
As a consequence, unlike K08, we detect 
no significant variability of the F606W 
flux between 2004 and 2006. Consider- 
ing the reduced and possible lack of vari- 
ability at F606W and the detection at 
F435W, several dust cloud models dis- 
cussed by K08 cannot be ruled out any- 
more (§ 14.3. 2p . K08 propose also a model of 
a Jovian planet surrou n ded b y a large disk 
of dust. IJanson et al.l ( 120121 ) exclude this 
explanation mainly because of the variabil- 
ity at F606W and the assumed dust belt 
crossing trajectory. Given our new pho- 
tometry and astrometry, we cannot reject 
this model (S l4T3T3j) . 

In the second part of our analysis, we 
study the possibility that Fomalhaut b 
is spatially resolved in our images. The 
signal-to-noise ratios of the detections are 
low and more data are required to con- 
firm the result but we find that our im- 
ages are more consistent with an extended 
source with diameter 1.16 AU than with a 
point source (§ 14.4.11 and 14.4.2"]) . The pho- 
tometric variability of an extended source 
model at F606W is larger than for a point 
source but it is not yet significant (< 2.5a, 
§ !4.4.3p . Two models are considered to ex- 
plain the size and the photometry of an 
extended source. First, the measurements 
are consistent with a cloud of dust pro- 
duced by a collision of two Kuiper belt ob- 
jects with radius 50 km that would have 
occurred ~ 50 years ago (§ 14.4.41) . The 
second model is an adaptation of the cir- 
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cumplane tary satellite swarm mo del pro- 
posed by iKennedyfeWyattl (j201ll ). It is 
consistent with the data when consid- 
ering a 2MEarth-lAf/ planet surrounded 
by a swarm that lies at 0.17-0.5 Hill 
radii f? ITO| . 

The nature of the Fomalhaut b ob- 
ject is still uncertain. However, from the 
two independent analysis of the HST data 
by K08 and our team, we can claim that 
Fomalhaut b is a real object that orbits 
Fomalhaut A. 
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